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Abstract 

Recently a cosmic ray propagation model has been introduced, where anisotropic diffusion is used as a 
mechanism to allow for O(IOO) km/s galactic winds. This model predicts a reduced antiproton background flux, 
suggesting an excess is being observed. We implement this model in GALPROP v50.1 and perform a analysis 
for B/C, ^"Be/^Be, and the recent PAMELA p/p datasets. By introducing a power-index parameter a that 
dictates the dependence of the diffusion coefficient D^x on height \z\ away from the galactic plane, we confirm 
that isotropic diffusion models with a — cannot accommodate high velocity convective winds suggested by 
ROSAT, while models with a — 1 (D^x oc |2:|) can give a very good fit. A fit to B/C and ^"Be/^Be data predicts 
a lower p/p fiux ratio than the PAMELA measurement at energies between approximately 2 GeV to 20 GeV. 
A combined fit including in addition the p/p data is marginal, suggesting only a partial contribution to the 
measured antiproton flux. 



1 Introduction 

In 2008-2009, the cosmic ray (CR) observation experiments PAMELA [1] and FERMI 
[2] both recorded an anomalous positron flux in excess of the expected astrophysical 
background. Similar excesses were also reported by the ATIC [3] and HESS [1] exper- 
iments at different energies. These observations initiated a flurry of activity aimed at 
explaining the discrepancy. Proposed solutions include both astrophysical sources, such 
as local pulsars O O [7] , as well as more exotic sources such as annihilating or decaying 
dark matter particles [H |9l [10] . The amount to which the latter contributes is extremely 
important information, as many models of new physics that extend the Standard Model 
predict the existence of particles that can make excellent candidates for dark matter. 

Inherent to any of the possible solutions is an assumption for the underlying physics 
describing how the CR particles are transported from their sources to Earth. At present, 
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our understanding of how CR propagate throughout the galaxy is mature but still in- 
complete. The most realistic description of CR propagation is currently obtained from 
models that assume CR transport occurs via a combination of spatial diffusion (resulting 
from interaction with the random and complicated galactic magnetic field) and convec- 
tion (resulting from interaction with large-scale streaming of CR away from the galactic 
disk). Explicit implementations have been presented in the literature. These include 
both semi-analytic [HI [T2] as well as fully numerical treatments [I3l Ull [151 IB] • 

In the canonical propagation scenario, the diffusion behavior is assumed to be identical 
everywhere in the galaxy: both in the gaseous disk as well as the surrounding region 
known as the halo. This isotropic diffusion model is capable of reproducing numerous 
astrophysical observations well, and this has led to its adoption as the model from which 
astrophysical backgrounds are normally derived. However, in isotropic diffusion models 
the convective wind velocity is limited to no more than about 10-20 km/s. Assuming a 
higher wind velocity results in the incorrect prediction of several existing observations, 
most notably measurements of secondary-to-primary ratios such as the boron-to-carbon 
(B/C) ratio, as well as the ratio of radioactive beryllium isotopes (^°Be/^Be) [T7 1 [T5 | [T^. 
These ratios, in particular, are highly sensitive to the parameters that control the various 
CR transport mechanisms, and are often used to benchmark predictions of a given model. 
While very high velocity galactic winds on the order of 1000 km/s have been observed 
in other galaxies [20], until recently it was thought that the Milky Way did not exhibit 
a wind, or that the wind velocity was limited to very low values. However, in 2007 the 
ROSAT satellite observed x-ray emission from within the Milky Way that is consistent 
with the presence of a galactic wind having a velocity of a few hundred km/s (2T1 122] . 
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The measurements are also consistent with a wind velocity profile that exhibits spatial 
dependence, and follows the radial distribution of supernova remnants (SNR) in the 
galaxy. The presence of high-velocity, radially dependent galactic winds appear to be 
natural and may also explain the large bulge-to-disk ratio observed by INTEGRAL 
[231121]. 

Recently, a propagation model has been introduced that employs both a radially de- 
pendent convective wind profile and anisotropic diffusion, where the diffusion behavior 
is not uniform throughout the galaxy, but instead varies as a function of position |25j . 
Specifically, the model assumes that the intensity of CR diffusion increases with distance 
from the galactic plane. This model is capable of supporting a convective wind veloci- 
ties of several hundred km/s, which is consistent with the ROSAT measurement, while 
simultaneously reproducing benchmark measurements such as the B/C and ^°Be/^Be 
ratios with similar accuracy to the traditional isotropic diffusion model. A crucial pre- 
diction of this scenario is that the astrophysical anti-proton background flux appears 
lower than what is predicted by isotropic diffusion models [26]. Such an outcome has 
significant implications for the indirect detection of particle dark matter, as it would 
reduce constraints on models for new physics that tend to over-predict the antiproton 
flux, assuming no excess is currently observed in the data [27] . 

In this paper we implement the anisotropic propagation model of [25] in the public CR 
propagation code GALPROP v50.1 [T3|, and explore varying levels of diffusion anisotropy 
by adjusting the rate at which the diffusion intensity increases with distance from the 
galactic plane. We perform a analysis for B/C, ^°Be/^Be and the recent PAMELA p/p 
data and find the models that fit the data best all exhibit spatially-dependent, anisotropic 
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diffusion. Furtlier, we find tliese models predict a.p/p flux tliat is significantly lower than 
the observed flux between 2 GeV and 20 GeV. When we attempt a combined flt including 
the p/p data, we flnd an increase in Xmim by A^^ = 35.4 for the 23 PAMELA data points 
considered. 

2 Anisotropic Diffusion 

As CR propagate they scatter and diffuse via interaction with turbulent fluctuations in 
the galactic magnetic fleld. Models of CR propagation traditionally assume diffusion 
is isotropic throughout the propagation volume. This implies a constant density of 
scattering centers throughout the volume, and an abrupt transition to free space (no 
scattering) at the outer boundary. In this case a propagating CR experiences, for a given 
momentum, the same scattering mean free path in the gaseous disk as it does far into 
the halo. In the presence of a galactic wind flow driving CR away from the gaseous 
disk and into the halo, the probability for a CR particle to return to the disk drops 
rapidly with increasing distance into the halo, as it not only must survive the convective 
action of the wind, but also an increasing number of scattering events on its return trip. 
The interplay between diffusive scattering and convective transport strongly affects the 
times a CR particle spends in the gaseous disk and halo regions of the galaxy. This, in 
turn, affects the prediction for the CR flux at Earth. In isotropic diffusion models the 
convective wind speed is limited to less than approximately 20 km/s before it is no longer 
possible to reproduce key observations such as the B/C ratio [TTl [13 [19] . This velocity 
limit is in conflict with measurements from the ROSAT satellite, which are consistent 
with galactic wind speeds up to several hundred km/s [211 [22] • 
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In the anisotropic propagation model of f25], CR diffusion is taken to be spatially 
dependent, and increases monotonically with distance above the gaseous disk. This 
corresponds to a gradual reduction in the density of scattering centers (and a smooth 
transition to free space) as CR travel from the disk towards the outer boundary of 
the propagation volume. In this scenario, the mean scattering length increases with 
increasing distance into the halo. Therefore, a CR particle in the halo that scatters back 
towards the disk can travel further before experiencing its next scattering event, which is 
likely in a region with a higher density of scattering centers. This "trapping" mechanism 
induces a drift of CR particles from of high diffusion to with low diffusion, 

which assists in counteracting the outward convection of CR due to the galactic wind. 
In this way, convection velocities of several hundred km/s can be accommodated while 
simultaneously preserving the CR residence times in both the disk and halo required to 
reproduce the CR fluxes observed at Earth. 

3 Cosmic Ray Transport Model 

For our study, we utilize the CR propagation code GALPROP v50.1 fTl]. GALPROP 
assumes that cosmic ray propagation in our galaxy is governed by the following transport 
equation 



where r,t) is the density of CR with total momentum p at location r at time t. 
The cosmic ray source distribution, as well as contributions from nuclear spallation and 
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decays are included in q{p,r,t). On the RHS, D^x is the spatial diffusion coefficient 
while V{r) is the velocity profile of an assumed galactic wind. The term proportional 
to Dpp accounts for random, stochastic accelerations (reacceleration) that result from 
scattering, and is equivalent to diffusion in momentum space. The term p(V ■ V)/3 
accounts for adiabatic loss or gain in momentum, in the case that the wind velocity is 
location dependent. Finally, the term with p = dp/dt accounts for momentum losses such 
as inverse Compton scattering, bremsstrahlung, as well as synchrotron radiation, while 
terms with r^^ and account for losses due to fragmentation and radioactive decay, 
respectively. GALPROP assumes a cylindrical propagation region having maximum 
radius from the galactic center r = rmax, and maximum extension above and below the 
galactic disk z = L. Free escape (\E' = 0) is assumed at the boundaries of the cylinder. 

3.1 GALPROP Modifications 

In this work the original GALPROP source code is modified in order to simulate anisotropic 
diffusion and a radially-dependent, ROSAT compatible convective wind. We follow 
closely the setup outlined in [25]. GALPROP utilizes the Crank- Nicholson implicit 
method [33] to numerically solve Eq. ([T]) and it is necessary to compute the finite- 
difference expansion coefficients in order to enact any modifications. The coefficients 
resulting from the inclusion of spatially dependent diffusion and convection have been 
computed and they agree to those in [25] in the limit of an equidistant computational 
grid. We assume the following grid spacing: Ar = 1 kpc and Az = 0.1 kpc through- 
out this analysis. The propagation cylinder height L is allowed to vary, however the 
maximum radius is fixed at rmax = 20 kpc. 

The spatial diffusion coefficient D^x is taken to be isotropic within a slab of thickness 
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1 kpc above and below the galactic plane, and anisotropic at distances > 1 kpc: 
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, \z\ < I kpc, (2) 

, l-^l > 1 kpc, (3) 

where (3 = v/c, Dq is a constant, p is the particle rigidity (momentum/electric charge), 
and 6 is the associated power index. The reference rigidity po is set to 4 GV. We introduce 
a power index a which governs the diffusion gradient; a = 1 was assumed in [25] and 
a = reproduces the original GALPROP parameterization. The diffusion coefficient is 
related to the reacceleration coefficient Dpp via 

D D =- ^ (4) 

where va is the speed of propagating Alfven waves: weak disturbances in the magnetic 
field that contribute to turbulent scattering. For the CR source distribution we adopt 
the parameterization of Case and Bhattacharya [34j: 

r — To \z\ 
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where vq = 8.5 kpc is the galacto-centric distance to Earth, and Zg = 200 pc. The distri- 
bution has a peak at r ^ 0.51ro. The galactic wind velocity is taken to be proportional 
to the source distribution at z = 

\V{r,z)\ =Q{r,z = 0){Vo + \z\-dV/dz). (6) 

We fix the constant velocity component at Vq = 100 km/s, and the gradient at dV/dz = 
35 km/s/kpc, which models the galactic wind velocity profile indicated by the ROSAT 
data. In GALPROP the galactic wind flow is one-dimensional along ±z and always 
outward from the galactic plane. The cosmic ray injection spectrum follows a power law 



in momentum, p~^. For nuclei, rj = 1.8 for rigidity p < 9 GV, and r] = 2.4 for p > 9 GV. 
Similarly, for electrons 7] = 1.6 for p < 4 GV, 2.54 for 4 GV < p < 100 GV, and 5.0 for 
p > 100 GV. We assume a constant Xqo factor (conversion factor from CO integrated 
temperature to H2 column density) of 1.9 x 10^°, and an H/He ratio of 0.11. 

4 Analysis 

Using these settings as a kernel, we scan parameters a, Dq, va, S, and L on the following 
grid: < a < 1 in steps of 0.2, 4 cm^/s < Dq < 6 cm^/s in steps of 0.2 cm^/s, 
30 km/s < Va < 60 km/s in steps of 2 km/s, 0.2 < 6 < 0.5 in steps of 0.02, and 
4 kpc < L < 8 kpc in steps of 1 kpc. For each model point {a, Dq, va, S, L} of the scan 
we compute the overall fo^' the resulting B/C, ^"Be/^Be, and p/p flux ratio curves. 
For our study we assume 

j=i i=i «i 

where index j runs over all M experimental datasets, each containing Nj measurements 
fij ± aij (i = 1 to Nj), and is the computed flux evaluated at the specific energy of 
the measurement i. Flux ratios are extracted using the CPLOT routine available on 
the GALPROP websitcl. For comparison with data, the local interstellar (LIS) flux was 
corrected for solar modulation using the "force field" approximation [35]. We assume 
the following modulation potentials, for B/C, ^'^Be/^Be, and p/p flux ratios, respectively, 
$ = 300 MV, 350 MV, and 600 MV. We include the following data. For B/C: HEAO-3 
[Sg, CREAM [37J, ACE [38j, and ATIC-2 [39]. For ^^Be/^Be: ISOMAX [IQ], and ACE 
For p/p: PAMELA 2010 [12]. Except for the PAMELA p/p measurements, all the 



^http: / /galprop. stanford.edu/web_galprop/galprop_homc. html 



data are obtained from the Galactic Cosmic Ray Database [13]. In 4 of the 23 PAMELA 
data points the reported uncertainty is asymmetric. For those points the larger error 
bar is used. No correlation among errors have been taken into account, as indicated by 
expression ([7]). Our exercise should therefore be considered as a first exploratory study of 
the allowed range of the parameters of the anisotropic diffusion model of CR propagation 
in our galaxy. 

4.1 Combined fit to B/C and ^^Be/^Be 

We first attempt a combined fit to both the B/C and ^°Be/^Be datasets. We find 
Xmin/'^dof = 38.3/31 within the range of the parameter space explored. The best-fit 
values and corresponding 3a bounds for parameters a, Dq, 6, and va are given in Table 
[11 A bound for L is not provided because the fit does not show significant dependence 
on the propagation cylinder height. 

Within the < a < 1 range of our scan, the Xmm occurs at the a = 1 boundary. 
The corresponding 3a allowed range is a > 0.27. The isotropic diffusion limit a = 
is excluded at 4cr, which confirms the finding of [25]. We remark that while our study 
focuses on the < a < 1 region of parameter space, a preliminary scan using a very 
coarse grid indicates that increases with a for a > 1. A detailed examination of the 
a > 1 region is left for a future analysis. 

The 3a allowed range of the diffusion constant Dq is 4.3 x 10^* cm^/s < -Do ^ 5.1 x 
10^^ cm^/s, while that of the diffusion index S is 0.35 < S < 0.43. These values are 
consistent with previous studies involving GALPROP [111[T3] (see also [T7]) and also with 
magneto-hydro dynamic theory, which predicts 6 = 1/3 for Kolmogorov-type turbulence 
power spectrum, and 6 = 1/2 for a Kraichnan-type spectrum. Similarly, the 3a range 
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for the Alfven velocity va is 47 km/s < f a < 58 km/s. These values are somewhat 
high compared to va — 30 km/s obtained for models with no convective wind pL3,J, and 
Va — 20 km/s expected for the warm ionized phase of the interstellar medium |45j . 
This increase is due to the large velocity gradient, dV/dz = 35 km/s/kpc, assumed for 
the galactic wind ([6]). A positive velocity gradient implies an adiabatic energy loss in 
response to the reduction in pressure of the CR fluid (term proportional to j9( V ■V)/3 in 
Eq. ([T])). The result is a redistribution of the B/C flux spectra towards lower energies. In 
order to reproduce the distinct peak in B/C at ~ 1 GeV, this reduction is compensated 
for with increased reacceleration, which is a mechanism that shifts CR from lower energy 
to higher energy. From Eq. (j4]) this translates into larger Dpp, or equivalently magnetic 
field disturbances that propagate with larger Alfven velocity (va)- 

Figure [1] shows the two-parameter likelihood surfaces for transport parameters Dq, 6, 
Va, and L as a function of a. Contours for Ax^ = 2.3, 4.61, 5.99, and 9.21 are given, 
corresponding to the 68%, 90%, 95% and 99% confidence level regions, respectively. 
The best-fit model point is indicated with a red diamond. Compared to models with 
isotropic diffusion (a = 0), those with spatially dependent, anisotropic diffusion (a > 0) 
provide significantly better fits to the data. The constraint on the propagation cylinder 
height L (Fig. [TJl) degrades with increasing a. As the diffusion gradient steepens, the 
distance at which a CR transitions into what is effectively free space (no scattering) 
occurs increasingly closer to the galactic plane, and away from the assumed boundary of 
the propagation region at L. This reduces the sensitivity of the model to the precise value 
of L, which no longer defines the extent of the magnetic corona, but becomes instead a 
computational parameter. In isotropic models, L strongly controls the residence time of 
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CR in the halo, as its position determines the location of the free escape boundary. 

Figure [2] shows the B/C and p/p flux ratio curves for model points that lie within 
the 68 % CL region, assuming a four parameter fit (Ax^ = 4.72). The grey (dark) lines 
correspond to the local interstellar (LIS) flux, while the green (light) lines correspond to 
the flux after correcting for solar modulation. The B/C flux is reproduced well. However, 
the p/p flux is significantly lower than the PAMELA measurement at energies between 
approximately 2 GeV and 20 GeV. 

4.2 Combined fit to B/C, ^^Be/^Be, and p/p 

We next attempt a combined fit to the B/C, ^°Be/^Be, and also the p/p data. We observe 
an increase in Xmin^ by Ax^ = 35.4 for the 23 PAMELA data points considered. Figure [3] 
shows the 68% CL curves obtained for the combined fit in red (dark) superimposed upon 
the results from Fig. [2] (B/C and ^''Be/^Be only) in green (light). The fit to the B/C 
data is now degraded somewhat as the models do not reproduce well the high precision 
HEAO-3 data above 1 GeV. The fit to the p/p data is marginally improved. However, 
an excess flux remains between approximately 2 GeV and 15 GeV. 

The lower half of Table [T] gives the best-fit parameter values and corresponding 3ct 
allowed range. The two-parameter likelihood surfaces are shown in Fig. |H Models with 
a > still provide the best fits to the data, however overall larger values of a (steeper 
diffusion gradient) are required: The best fit occurs at the a = 1 boundary of our scan 
range (0 < a < 1), and the 3a allowed range is a > 0.51. The allowed range for 
parameters Dq, va, and S are now shifted to lower values relative to those obtained for 
the fit to B/C and ^^Be/^Be only, with the most pronounced reduction occurring for 
the Alfven velocity. A reduction in Alfven velocity corresponds to weakening of the 
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reacceleration mechanism (smaller Dpp in (jl])), a process that boosts CR momentum to 
higher values. Weaker reacceleration causes both the B/C and p/p flux distributions to 
shift to lower energies, which improves the fit to p/p but simultaneously degrades the fit 
to B/C. Lower values for Dq and 6 reduce the intensity of spatial diffusion, which forces 
CR to remain longer in the gaseous disk. This increases the rate of secondary production 
and boosts the B/C flux, which helps to compensate for the reduction in Alfven velocity. 

5 Conclusions 

We consider the anisotropic propagation model of [25], which assumes that CR diffusion 
increases linearly with distance from the galactic disk, and radially dependent convective 
wind flows with maximum velocity of O(IOO) km/s suggested by the ROSAT measure- 
ment. This model is capable of reproducing key observations such as the B/C and Be/Be 
flux ratios, whereas traditional models based on isotropic diffusion cannot support con- 
vective wind velocities of more than approximately 20 km/s. 

We implement this model in GALPROP, and introduce a new parameter, a, which 
controls the dependence of the diffusion coefficient D^^ on the distance \z\ above the 
gaseous disk; a = corresponding to isotropic diffusion, while a = 1 a linear increase as 
in [25]. We then perform a analysis for the B/C, ^°Be/^Be, and the recent PAMELA 
p/p data sets in order to constrain the transport parameters. For a combined fit to the 
B/C and ^°Be/^Be flux only, we observe Xmin/'^-dof = 38.3/31, indicating a good fit to both 
data sets. We find a > 0.27, 4.3 x lO^^ cmVs <Do<5.lx lO^^ cmVs, 0.35 <S< 0.43, 
and 47 km/s < fyi < 58 km/s at 3a for the range of parameter space we explored. 
The Alfven velocity is higher than expected, and due to a large wind velocity gradient, 
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dV/dz = 35 km/s/kpc, required to model the wind profile suggested by the ROSAT 
measurements. The isotropic diffusion case a = is excluded at Aa. Models within 
the 68% CL region predict a lower p/p flux than the PAMELA observation between 2 
GeV to 20 GeV. For a combined fit to B/C, ^^Be/^Be and also the PAMELA p/p data 
we observe an increase in Xmin' 

by Ax^ = 35.4 for the additional 23 PAMELA data 
points considered, indicating the overall fit is marginal. In this case we find a > 0.51, 
4.1 X 10^8 cmVs <Do< 4.8 x lO^^ cmVs, 0.33 <6 < 0.38, and 43 km/s < < 49 km/s 
at 3cr, which are consistent with expectations. Models within the 68% CL region give 
an improved fit to the p/p flux, however the predicted distributions are still below the 
PAMELA observation between approximately 2 GeV to 15 GeV. 

These results suggest that spatially dependent diffusion represents a viable solution 
to the problem of accommodating (9(100) km/s convective wind flows in models of CR 
propagation, and that in this case the observed antiproton flux may exhibit an excess 
above the astrophysical background. An anomalous excess in the antiproton flux would 
have many implications for the indirect detection of particle dark matter (DM). For ex- 
ample, pair annihilation of light neutralino DM can explain the positron excess observed 
by e.g. PAMELA ^6]. These annihilations necessarily produce hadrons in addition to 
leptons, and so a contribution to the antiproton flux is also expected from this scenario. 
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Combined Fit 


Parameter 


Best Fit Value 


3(7 Range 




a 


1.0 


> 0.27 


B/C, lOBe/^Be 


Dn 


4 fi X 10^^ rm^ /s 


4 3-51 xin^^ rm^ /s 




VA 


52 km/s 


47 - 58 km/s 




5 


0.38 


0.35 - 0.43 




a 


1.0 


> 0.51 


B/C, lOBe/^Be, p/p 




4.4 X 10^8 cm^ /s 


4.1 - 4.8 xlO^^ cm^/s 




VA 


46 km/s 


43 - 49 km/s 




S 


0.36 


0.33 - 0.38 



Table 1: Best-fit model parameters for the combined fits explored in this study. The 3cr allowed range is given 
for L = 5 kpc for the fit to B/C and ^^Be/^Be (above), and for L = 6 kpc for the fit including p/p (bottom), where 

we find the minimum x^- No bound is provided for parameter L, as the fit does not exhibit strong dependence on 
the propagation cylinder height. 
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(a) Diffusion Constant vs a 



o 

o 
Q 




(b) 5 vs a 



(c) Alfven Velocity vs a 





(d) Propagation Cylinder Height vs a 




Figure 1: Two-parameter likelihood surfaces for the transport parameters as a function of the diffusion gradient 
index a, obtained for combined fit to B/C and ^°Be/^Be flux ratios. Contours for Ax^ = 2.3, 4.61, 5.99, and 
9.21 are given, corresponding to the 68%, 90%, 95% and 99% confidence level regions, respectively. Subfigures: 
(a) diffusion constant Dq, (b) diffusion coefficient power index S, (c) Alfven velocity va, (d) propagation cylinder 
maximum height L. The red diamond indicates the location of xLin- Whereas isotropic diffusion models with 
a = {Dxx independent of position) cannot accommodate the high velocity convective winds suggested by ROSAT, 
models with a = 1 {Dxx oc \z\) can give a very good fit. 
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Figure 2: Predicted fiux distribution for models within the 68% confidence level region (four parameter fit 
Ax^ = 4.72), for the combined fit to B/C and ^''Be/^Be. The grey (dark) curves indicate the local interstel- 
lar (LIS) flux, while the green (light) curves have been corrected for solar modulation using the "force field" 
approximation (modulation potential indicated). The B/C fiux ratio is reproduced well, however the predicted p/p 

fiux is significantly lower than the PAMELA measurements between approximately 2 GeV to 20 GeV. 
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Figure 3: The predicted flux distribution obtained for models within the 68% confidence level region for the combined 

fit to the B/C, -'^°Be/^Be, and p/p data in red (dark) superimposed on the results obtained for the fit to B/C and 

^"Be/^Be only in green (light). The xL\n observed to increase by Ax^ = 35.4 for the additional 23 PAMELA 

data points considered. The fit to the B/C flux is degraded as the distribution does not reproduce well the high 

precision HEAO-3 measurements above 10 GeV. The fit to p/p is marginally improved, however the distribution is 

20 

still lower than the PAMELA observation between approximately 2 GeV and 15 GeV. 



(a) Diffusion Constant vs a 



o 

Q 




(b) 5 vs a 



(c) Alfven Velocity vs a 





(d) Propagation Cylinder Height vs a 




Figure 4: Two-parameter likelihood surfaces for the transport parameters as a function of the diffusion gradient 
index a, obtained for combined fit to B/C, ^°Be/^Be and p/p flux ratios. The contours are as in Fig. [TJ We find 
an increase of Xmin = 35-4 for the 23 PAMELA data points considered, and a significant reduction of the 

Alfven wave velocity {va) 
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